Matrix valued Brownian motion and a paper by Polya 

Philippe Biane 

1. Introduction 

This paper has two parts which are largely independent. In the first one I 
recall some known facts on matrix valued Brownian motion, which are not so easily 
found in this form in the literature. I will study three types of matrices, namely 
Hermitian matrices, complex invertible matrices, and unitary matrices, and try to 
give a precise description of the motion of eigenvalues (or singular values) in each 
case. In the second part, I give a new look at an old paper of G. Polya |14j, where 
he introduces a function close to Riemann's ^ function, and shows that it satisfies 
Riemann's hypothesis. As put by Marc Kac in his comments on Polya's paper [llj . 
"Although this beautiful paper takes you within a hair's breadth of Riemann's 
hypothesis it does not seem to have inspired much further work and reference to 
it in the mathematical litterature are rather scant". My aim is to point out that 
the function considered by Polya is related in a more subtle way to Riemann's 
^ function than it looks at first sight. Furthermore the nature of this relation is 
probabilistic, since these functions have a natural interpretation involving Mellin 
transforms of first passage times for diffusions. By studying infinite divisibility 
properties of the distributions of these first passage times, we will see that they 
are generalized gamma convolutions, whose mixing measures are related to the 
considerations in the first part of this note. 

2. Matrix Brownian motions 

We will study three types of matrix spaces, and in each of these spaces consider 
a natural Brownian motion, and show that the motion of eigenvalues (or singular 
values) of this Brownian motion has a simple geometric description, using Doob's 
transform. The following results admit analogues in more general complex sym- 
metric spaces, but for the sake of simplicity, discussion will be restricted to type 
A symmetric spaces. Actually the interesting case for us in the second part will 
be the simplest one, of rank one, but I think that this almost trivial case is better 
understood by putting it in the more general context. Some references for results 
in this section are [2], |4], [5], [7], [8], [9], [10], [13], [H], [16] . 

2.1. Hermitian matrices. Consider the space of n x n Hermitian matrices, 
with zero trace, endowed with the quadratic form 



{A,B) = Tr{AB). 
1 
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Let M(t) be a Brownian motion with values in this space, which is simply a Gaussian 
process with covariance 

E[Tr{AM{t))Tr{BM{s))] = Tr{AB)s A t 

for A, B traceless Hermitian matrices. 

Let Xiit) > X2{t) > ... > Xn{t) be the eigenvalues of M{t); they perform a 
stochastic process with values in the Weyl chamber 

e = {{xi,...,Xn) G R" Ixi > a;2 > ... > a;„}nF„ 

where 

n 

= I (a;i, . . . , a;„) G R" | ^ = | . 

i=l 

Let P( be the transition probability semi-group of Brownian motion killed at the 
boundary of the cone C. This cone is a fundamental domain for the action of the 
symmetric group Sm which acts by permutation of coordinates on R". Using the 
reflexion principle, one shows easily that 

where Pt{x,y) = (27ri)~*^"~-^'/^e~l^~^l^/^* and e(c7) is the signature of a. Let h be 
the function 

i>j 

Proposition 2.1. The function h is the unique (up to a positive multiplicative 
constant) positive harmonic function for the semigroup p^, on the cone C, which 
vanishes on the boundary. 

The harmonic function h corresponds to the unique point at infinity in the 
Martin compactification of C. Consider now the Doob's transform of which is 
the semigroup given by 

It is a diffusion semigroup on C with infinitesimal generator 

iA + (Vlog/i,V-) . 

Proposition 2.2. The eigenvalue process of a traceless Hermitian Brownian 
motion is a Markov diffusion process in the cone C, with semigroup g*. 

We can summarize the last proposition by saying that the eigenvalue process is 
a Brownian motion in C, conditioned (in Doob's sense) to exit the cone at infinity. 

2.2. The group 5L„(C). This is the group of complex invertible matrices 
of size n X n, with determinant 1. Its Lie algebra is the space s[„(C) of complex 
traceless matrices. Consider the Hermitian form 

{A,B) = Tr{AB*) 

on s[„(C) which is invariant by left and right action of the unitary subgroup SU (n). 
This Hermitian form determines a unique Brownian motion with values in sln{C). 
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The Brownian motion gt, on S'L„(C), is the stochastic exponential of this Brownian 
motion, solution to the Stratonovich stochastic differential equation 

dgt = gtdwt 

where Wt is a Brownian motion in 5[„(C). 

There are two remarkable decompositions of S'L„(C), the Iwasawa and Cartan 
decompositions. The first one is 5'L„(C) = NAK where K is the compact group 
SU{n), A is the group of diagonal matrices with positive coefficients, and determi- 
nant one, and N is the nilpotent group of upper triangular matrices with I's on 
the diagonal. Each matrix of SLn{C) has a unique decomposition as a product 
g — nak of elements of the three subgroups N,A,K. This can be easily inferred 
from the Gram-Schmidt orthogonalization process. If gt is a Brownian motion in 
SLn{C), one can consider its components nt,at,kt. In particular, denoting by 
(e"'^'^*-', . . . , e"""*^*-*) the diagonal components of at the following holds (cf TS]). 

Proposition 2.3. The process . . . , w„(i)) is a Brownian motion with 

a drift p = {—n + 1, — n -I- 3, . . . , n — 1) in the subspace iJ„. 

The other decomposition is the Cartan decomposition SLn{C) — KA^K, 
where A'^ is the part of A consisting of matrices with positive nonincreasing coef- 
ficients along the diagonal. In order to get the Cartan decomposition of a matrix 
g € SLn{C), take its polar decomposition g = ru with r positive Hermitian, and 
u unitary, then diagonalize r which yields g = vav' with v and v' unitary and a 
diagonal, with positive real coefficients which can be put in nonincreasing order 
along the diagonal. These coefficients are the singular values of the matrix g. This 
decomposition is not unique since the diagonal subgroup of SU (n) commutes with 
A, but the singular values are uniquely defined. Call (e^^^^'^ . . . , e""*^*-') the singular 
values of the Brownian motion gt, with oi > 02 > . . . > a„. They form a process 
with values in the cone C. Let us mention that this stochastic process can also be 
interpreted as the radial part of a Brownian motion with values in the symmetric 
space S Ln{C) / SU (n) . We will now give for the motion of singular values a similar 
description as the one of eigenvalues of the Hermitian Brownian motion. For this, 
consider a Brownian motion in i?„, with drift p, killed at the exit of the cone C. 
This process has a semigroup given by 

p°-''{x,y)^e'''''y-^^-'^'''''^/Vt{x,y). 

Proposition 2.4. The function 

h^iy) = [1(1 -e^fe-J'')) 

is a positive harmonic function for the semigroup p}''^ , in the cone C, and vanishes 
at the boundary of the cone. 

It is not true that this function is the unique positive harmonic function on the 
cone; indeed the Martin boundary at infinity is now much larger and contains a 
point for each direction inside the cone, see The Doob-transformed semigroup 

is a Markov diffusion semigroup in the cone C, with infinitesimal generator 

^A+{p,-) + {Vloghfy-). 
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Note that it can also be expressed as 

iA+(Vlog/^^V•) 

with 

h^iy) ^Yisinhiyj ~ y^) 

i>j 

(see [10]). 

Proposition 2.5. The logarithms of the singular values of a Brownian motion 
in SLn{C) perform a diffusion process in the cone C, with semigroup . 

As in the preceding case, we can summarize by saying that the process of 
singular values is a Brownian motion with drift p in the cone C, conditioned (in 
Doob's sense) to exit the cone at infinity, in the direction p. 

2.3. Unitary matrices. The Brownian motion with values in SU{n) is ob- 
tained by taking the stochastic exponential of a Brownian motion in the Lie algebra 
of traceless anti-Hermitian matrices, endowed with the Hermitian form 

{A,B) = -Tr{AB). 

Let e*^% . . . , e*^" be the eigenvalues of a matrix in SU{n), which can be chosen so 
that = 0; s-i^d Oi > 92 > ■ ■ ■ > On, di — On < Stt. These conditions determine 

a simplex A„ in which is a fundamental domain for the action of the affine 
Weyl group on i/„. Recall that the affine Weyl group W is the semidirect product 
of the symmetric group Sm which acts by permutation of coordinates in i7„, and 
of the group of translations by elements of the lattice (27rZ)" n Hn- 

One can use the reflexion principle again to compute the semigroup of Brownian 
moton in this simplex killed at the boundary. One gets an alternating sum over the 
elements of 

wew 

The infinitesimal generator is 1/2 x the Laplacian in the simplex, with Dirichlet 
boundary conditions. It is well known that this operator has a compact resolvent, 
and its eigenvalue with smallest module is simple, with an eigenfunction which can 
be chosen positive. Consider the function 

Proposition 2.6. The function /i" is positive inside the simplex A„, it van- 
ishes on the boundary, and it is the eigenfunction corresponding to the Dirichlet 
eigenvalue with smallest module on A„. This eigenvalue is A = {n — n^)/6. 

The Doob-transformcd semigroup 

is a Markov diffusion semigroup in A„, with infinitesimal generator 

iA + (Vlog;i",V-)-A. 



MATRIX VALUED BROWNIAN MOTION AND A PAPER BY p6lYA 



5 



Proposition 2.7. The process of eigenvalues of a unitary Brownian motion is 
a diffusion with values in A„ with probability transition semigroup q^. 

Again a good summary of this situation is that the motion of eigenvalues is 
that of a Brownian motion in the simplex A„ conditioned to stay forever in this 
simplex. 

2.4. The case of rank 1. In the next section we will need the simplest case, 
that of 2 X 2 matrices. Consider first the case of Hermitian matrices. The pro- 
cess of eigenvalues is essentially a Bessel process of dimension 3, with infinitesimal 
generator 

I (f I d 

obtained from Brownian motion killed at zero, of infinitesimal generator 

with Dirichlet boundary condition at 0, by a Doob transform with the positive 
harmonic function h{x) = x 

In the case of the group 5*^2(0), or the symmetric space SL2{C)/SU{2), which 
is the hyperbolic space of dimension 3, the radial process has infinitesimal generator 

1 d 

obtained from Brownian motion with a drift 

1 ^ d 

2 dx'^ dx 

with Dirichlet boundary condition at 0, by a Doob's transform with the function 
1 - 

Finally the last case is Brownian motion in SU{2), where the eigenvalue process 
takes values in [0, n] and has an infinitesimal generator 

1 rf2 ^^^^ d 

2 dx'^ dx 

obtained by a Doob transform at the bottom of the spectrum from 

2'dx^ 

on [0,7r] with Dirichlet boundary conditions at and tt, by the function sin(a;). 

For these last two examples, we shall write a spectral decomposition of the 
generator Lj, i = 1, 2, of the form 

(2.1) f{x) = I ^\{x) [ I ^\{y)f{y) drmiy)] di^iW i = 1,2 

for every / e L'^{mi), where mj is measure for which Li is selfadjoint in L^{mi), 
and the functions $\ solutions to 

Li^{ + = 

and Vi is a spectral measure for Li. 

For Li = on [0, tt] with Dirichlet boundary conditions, mi (dec) is Lebesgue 
measure on [0, tt], and Li is selfadjoint on L^([0,7r)]. Furthermore 

^\{x) = sm{^/2\x) 
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and 

oo 

(2.2) MdX) = -J2^ny2{dX). 

n=l 

For L2 — ^-j^ + [0; +00 [; the measure m2{dx) = e^^dx, and 

$i(a;) = e-=" sin(\/2A - Ix) A > 1/2. 
The spectral measure is 

(2.3) MdX)^ — dA A>l/2 

TTv 2A — 1 

on [l,+oo[. Of course formulas (|2.ip . (|2.2|) . (|2.3p are immediate consequences of 
ordinary Fourier analysis. 

Note that the spectral decompositions, and in particular the measures de- 
pend on the normalisation of the fonctions <f>A- We have made a natural choice, 
but it does not coincide with the usual normalisation of Weyl-Titchmarsh-Kodaira 
theory, see [6]. 



3. MacDonald's function and Riemann's ^ function 

3.1. Polya's paper. In his paper |14| . Polya starts from Riemann's ^ function 

f(s) = s(s-l)7r-^/2r(s/2)C(s) 

where ^ is Riemann's zeta function. Then ^ is an entire function whose zeros are 
exactly the nontrivial zeros of ^. 
Putting s = 1/2 + jz yields 

/>oo 

^(z) = 2 / $(u) cos{zu)du 
Jo 

with 

00 

(3.1) $(w) = 27re^"/2 ^{2TTe^''n^ - 3)n2e-™''='" 

n=l 

and the function $ is even, as follows from the functional equation for Jacobi 9 
function; furthermore 

so that 

$(u) ^ 4^2(g9n/2 _^ g-9«/2)g-.(e-+e--) ^ ^ 

This lead Polya to define a "falsified" ^ function 

C (z) = 87r2 / (e9"/2 + g-9«/2)g-.(e-"+e--) cos(zu)d7.. 

Jo 

The main result of 14| is 

Theoreme 3.1. The function ^* is entire, its zeros are real and simple. Let 
N{r), (resp. N*{r)) denote the number of zeros of ^{z) (resp. £,*{z)) with real part 
in the interval [0,r], then N{r) — N*(r) = O(logr). 
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Recall that the same assertion about the zeros of the function ^ (without the 
statement about simplicity, beware also that s = 1/2 + iz) is Riemann's hypothesis. 
Recall also the well known estimate 

7V(r) = f log(r/27r)-f + 0(1). 
zvr zvr 

Polya's results rely on the intermediate study of the function 

/■OO 



from which ^* is obtained by 

Ciz) = 27r2(©(iz/2 - 9/4, tt) + (5(iz/2 + 9/4, tt)) 

Polya shows that &{z, a) has only purely imaginary zeros, (as a function of z) and 
the number of these zeros with imaginary part in [0, r] grows as ^ log ^ — + 0{1). 
The results on ^* are then deduced through a nice lemma which played a role 
in the history of statistical mechanics (the Lee- Yang theorem on Ising model) , as 
revealed by M. Kac [llj . We shall now concentrate on &(z,a). In particular, for 
a = TT, the function ^(z) = (5(iz/2,7r) is another approximation of ^ which has 
many interesting structural properties. 

3.2. MacDonald functions. The function denoted 6(z, a) by Polya is actu- 
ally a Bessel function. Indeed, MacDonald's function, also called modified Bessel 
function (see e.g. [1]), given by 



/•OO 

Jo 



satisfies Kz{2x) = (&{z,x). The function (&{z,a) is therefore essentially a MacDon- 
ald function, as noted by Polya. MacDonald function is an even function of /i and 
satisfies 

(3.2) —K,,{x) = K^+i{x) - K^^iix) 

X 



(3.3) - 2—K^{x) = K^+,(x) + K^_,{x) 



d 

dx 

The first of these equations is used by Polya in a very clever way to prove that the 
zeros (in z) of 25(z, a;) are purely imaginary. 



3.3. Spectral interpretation of the zeros. From 
from which one gets 

= (_i£±l _ d _ d 

^ X dx ' ^ X ax 

= (JL + 1 a _ !L.)K 

This differential equation will give us a spectral interpretation of the zeros of ©(z, x) 
Change variable by ipf_i{x) — K^{e^) to get 

(3.4) (_^ + e2-)^^ = -/.V;. 
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Since Kf^ vanishes exponentially at infinity, the spectral theory of Sturm-Liouville 
operators on the half- line (see e.g. [6], [12j ) implies that the squares of the zeros 
oi fjb ^ ipfj,{y) are the eigenvalues of — e^^ on the interval [y, -|-oo[ with the 
Dirichlet boundary condition at y, the functions being the eigenfunctions. Since 
this operator is selfadjoint and negative the zeros are purely imaginary, and are 
simple. 

This spectral interpretation of the zeros of MacDonald function is well known 
[17j . I do not know why Polya does not mention it. 



3.4. H — xp. Equation (13. 4p can be put into Dirac's form, indeed the equa- 
tions 

(-£ + 1 + ^^1.9 -7/ 

imply 

Using the change of variables u = , we get 

+ =75 

Remark that this Dirac system yields a perturbation of the Hamiltonian H — xp 
considered by Berry et Keating [3], in relation with Rieniann's zeta function. 

3.5. Asymptotics of the zeros. General results on Sturm-Liouville opera- 
tors allow one to recover the asymptotic behaviour of the spectrum, thanks to a 
semiclassical analysis, see e.g. |12| . One can get a more precise result using the 
integral representation of Ki^. Polya gives the asymptotic estimate 



rZiry 



0(e-f!/yl^l-3/2) 



in the strip |a:| < 1 uniformly as y — > cxd, where 

$ = ylog- - y- - . 

a 4 

This estimate can be obtained by the stationary phase method, writing 

/•oo 



Making a contour deformation we get 

_j_ gZ(t-(-if )-2aisinh(t)^^ _ ^ J'^ /'^ ^z{A-it+i^)-2a cosh(A-it+i^) 

and Polya's estimate can be obtained by standard methods, which give also esti- 
mates for the derivatives of MacDonald's function. Finally the zeros of y ^ Kiy{2a) 
behave like the solutions to 

. y TT 1 

V log V = (n -\ — vr n mteger 
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The number of zeros with imaginary part in [0, T] is thus ^ log ^ — ^ + 0{1). 

4. Probabilistic interpretations 

We will now give interpretations of the functions £^ and using first passage 
times of diffusions. 

4.1. Brownian motion with a drift. The first passage time at a; > of 
Brownian motion started at follows a 1/2 stable distribution i.e., 



P(T^. e dt) ^ T. " d± 



with Laplace transform 



e 



Adding a drift a > to the Brownian motion gives a first passage distribution 
with Laplace transform 

This is a generalized inverse Gaussian distribution. In particular, its Mellin trans- 
form is 

K_i/2[ax) 

which gives a probabilistic interpretation of MacDonald's function (as a function of 
s) as a Mellin transform of a probability distribution. 

4.2. Three dimensional Bessel process. There exists a similar interpreta- 
tion of the ^ function, which is discussed in details in |4], [5], for example. Consider 
the first passage time at a > of a three dimensional Bessel process (i.e., the norm 
of a three dimensional Brownian motion) starting from 0. The Laplace transform 
of this hitting time is 

sinh Xa 

Let S'^ be an independent copy of Sa, and let 

Wa=Sa+S',; 

then the density of the distribution of Wa is obtained by inverting the Laplace 
transform. One gets 

oo 

P{Wa e dx) = ^{-K^n^x/a^ - 37r2nVa2)e~'^'"'^/2"'dx 

n=l 

from which one can compute the Mellin transform 

E[W:] = 2(2aV^)'f(2s). 
The function 2^ thus has a probabilistic interpretation, as Mellin transform of 
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4.3. Infinite divisibility. The distributions of Tx and Wa are infinitely di- 
visible. Indeed 

log£;°[exp(-^r^)] = -xTA^T^+aj; 
which shows that is a subordinator with Levy measure 



Similarly 

log£;[exp(-^l^a)] = 21og(Aa/sinh(Aa)) 
therefore the variable Wa has the distribution of a subordinator, with Levy measure 

oo 

E2 2 J. / 2 
g-^ n t/a 

n=l 

taken at time I. Observe however that the process {Wa)a>o is not a subordinator. 
4.4. Generalized gamma convolution. The gamma distributions are 

Ph^,c e dt) = ^r~^e-'/'dt = r^Adt) 
r{uj) 

where to and c are > parameters. The Laplace transform is 

^[e-A7.,e] = (1 + A/c)-". 

The gamma distributions form a convolution semigroup with respect to the param- 
eter uj, i.e., 

^U-L.C * L^2,C ^UJI-{-U2tC • 

The Levy exponent of the gamma semigroup is 

poo -ct 

V'c(A) = log(l + A/c) = / (1 - e-^*) — 
Jo 

so that this is the semigroup of a subordinator with Levy measure e~'^*/tdt. 

The generalized gamma convolutions are the distributions of linear combina- 
tions, with positive coefficients, of independent gamma variables, and their weak 
limits. 

One can also characterize the generalized gamma convolutions as the infinitely 
divisible distributions with a Levy exponent of the form 

/>oo 

V'(A) = / M^)dHc) 
Jo 

for some positive measure p which integrates I/c at oo. This measure is called the 
Thorin measure of the generalized gamma distribution. The variables and Wa 
of the preceding paragraph are generalized gamma convolutions. Indeed it is easy 
to check, using the computations of section l473l that Wa has a generalized gamma 
convolution as distribution, with Thorin measure 

oo 

(4.1) i^{dc) = 2Y,5nya-{dc) ; 

n=l 
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whereas Tx is distributed as a generalized gamma convolution with Thorin measure 

dc 



fir 

(4.2) iy{dc) = c > aV2 



since 



-aH/2 /.oo _ 

e 



a2/2 TT^C - 0^/2 

4.5. Final remarks. We can now make a connection between the preceding 
considerations and those of the first part of the paper. Indeed, the Thorin measures 
associated with the variables and Wa can be expressed as spectral measures 
associated with the generators of Brownian motion on matrix spaces. The hitting 
times of Brownian motion with drift are related with the radial part of Brownian 
motion in the symmetric space S L2{C) / SU (2) , whereas the hitting times of the 
Bessel three process are related with the Brownian motion on the unitary group 
SU{2). The precise relations are contained in formulas ((^ . ((O)) . ((i7T|) . (jO)) . 
Thus the Riemann ^ function, which is the Mellin transform of a hitting time of 
the Bessel three process, as in section 4.2, and the Polya ^ function from section 
3.1, which appears as Mellin transform of hitting time of Brownian motion with 
drift, are related in this non obvious way. 
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